INSTITUTE  FOR  TELECOMMUNICATION  SCIENCES  BOULDER  CO  F/S  20/14 

PROPAGATION  IN  A  TROPOSPHERIC  DUCT  WITH  A  SINGLE-STEP  DISCONTIN— ETC ( U) 
JUN  80  R  H  OTT  MlPR-FYl 17579N21 19 

AFWAL-TR-80-1138  NL 


UNCLASSIFIED 


1  OF  | 

■ 

■ 

■ 

■ 

■ 

| 

■ 

■ 

1 

■ 

■ 

■ 

■ 

END 

1 1-80 
DTIC 

A090544 


PROPAGATION  IN  A  TROPOSPHERIC  DUCT  WITH  A  SINGLE-STEP  DISCONTINUITY  IN  THE 
REFRACTIVE  INDEX  IN  THE  DIRECTION  OF  PROPAGATION 


R.  H.  Ott 

Institute  for  Telecommunication  Sciences 
Department  of  Commerce  j 

Boulder,  Colorado  80303 


June  1980 


DTIO 


r 

V 


—  .e:cte| 

OCT  1  7  1980 


1 

1 


TECHNICAL  REPORT  AFWAL-TR- 80-1138 


Final  Report  for  Period  Novemberv1979  -  May  1980 


Approved  for  Public  Release;  distribution  unlimited 


AVIONICS  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
MRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


80 


1  r\ 
j  u 


14  0 


NOTICE 


When  Government  drawings,  specifications ,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related  Government  procurement  operation 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  the  fact  that  the  government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  re¬ 
garded  by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 


other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture 
use,  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ASD/PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS,  it  will 
be  available  to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


RAYMOND  P.  WASKY,  Project  Engineer 
EO  Sensor  Science  &  Engineering  Group 
EO  and  Reconnaissance  Branch 
Avionics  Laboratory 


HAROLD  E.  GELTMACHER,  Chief 
EO  Sensor  Science  &  Engineering  Group 
Recon  &  Weapon  Delivery  Division 
Avionics  Laboratory 


FOR  THE  COMMANDER t 


Electro-Optics  Recon  Branch 
Reconnaissance  §  Weapon  Delivery  Division 
Avionics  Laboratory 


"If  you r  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list,  or 
if  the  addressee  is  no  longer  employed  by  your  organisation  please  notifu  avwat./aart, 
W-PAFB,  OH  45433  to  help  us  maintain  a  current  mailing  list". 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security 
considerations ,  contractual  obligations,  or  notice  on  a  specific  document. 

AIR  RORCC/M7M) /**  tap  tarn  Mr  1M0  -  »0 

I 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


ril.  REPORT 


(iff)  REPORT  DOCUMENTATION  PAGE 


JJ  ^PROPAGATION  IN  A  TROPOSPHERIC  JUJCT  WITH  A  TINGLE-1 
r  iTEP  DISCONTINUITY  IN  THE  REFRACTIVE  INDEX  IN  TH 

Direction  of  propagation . — - 


Ml 


Ty/uMBER 

■^TR-ajf)-n38~ 


VTITI  £  Cmrt  ZiAdtlm) 


MW  &A  & _ 


\  cnin  r  Fr-; 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


:atalog  number 


T VP F  OF  REPORT  >  PERIOD  COVERED 


7.  author^; 


0tt 


(/£  MIPR' 


/On- 


9.  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Institute  for  Telecommunication  Sciences 
Department  of  Commerce 
Boulder,  Colorado  80303 

II.  CONTROLLING  OFFICE  NAME  AND  AOORESS  (/  /J 

Air  Force  Avionics  Laboratory  (AAR)  v--'i 

Air  Force  Wright  Aeronautical  Laboratories,  AFSC 
Wright-Patterson  Air  Force  Base,  Ohio  45433 


1 4.  MONITORING  AGENCY  NAME  A  ADDRESS^//  different  from  Controlling  Office) 


Final  EeF&jt.ior  Period- 
79-May 


9.  CONTRACT  OR  GRANT  NUMBERfaJ 


FY1 1 7579N21 1 9 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 

- AJIOajMUilT  I - 


( ( (*\  Project  ILIR  I  Sl 

—  Task  ILIR52  ^ - 


Work  Unit  II.IR9209 


12.  REPORT  DATE 


13.  NUMBER  OF  PAGES 


26 


15.  SECURITY  CLASS,  (of  thla  report) 

UNCLASSIFIED 


15a.  DECLASSIFICATION  DOWNGRADING 
schedule 


16.  DISTRIBUTION  STATEMENT  (of  thla  Report) 

Approved  for  public  release;  distribution  unlimited. 


J  47.  DIS  -<IBUTION  STATEMENT  (of  the  abatrect  entered  In  Block  20,  If  different  from  Report) 


IS.  SUPPLEMENTARY  NOTES 


f £ ;  j  f  p 


19.  KEY  WORDS  (Confirm*  on  reveree  aide  II  neceaaary  and  Identity  by  block  number) 


20.  ABSTRACT  (Confirm*  on  revoroa  aide  If  neceaaary  and  Identify  by  block  number) 


A  Green's  function  approach  is  used  to  derive  an  expression  for  the  field 
inside  a  laterally  Inhomogeneous  duct.  The  laterally  inhomgeneous  duct  is 
assumed  to  have  a  single  step,  and  convergence  criteria  for  the  step  size 
and  number  of  modes  are  discussed. 


DC 

1  JAN  71  1473  COITION  OP  '  NOV  «•  I*  OOSOLETC 


1ICURITV  CLASSIFICATION  OP  THIS  PAGE  (Whin  Dm  gnu, id) 


Vo  3  *>vi 


A' 


r 


FOREWORD 


This  technical  report  describes  the  work  completed  by  the  Institute  for 
Telecommunication  Sciences,  U.S.  Department  of  Commerce,  Boulder,  Colorado,  under 
Project  ILIR9209,  "Over-the-Horizon  Target  Detection  Feasibility  Study."  This 
work  has  been  supported  by  the  Avionics  Laboratory  Director's  Fund. 

The  work  described  herein  is  for  the  period  November  1979  to  May  1980,  under 
the  direction  of  Mr.  Raymond  P.  Wasky  (AFWAL/AARI-3) ,  Electro-Optics  and  Recon¬ 
naissance  Branch,  Reconnaissance  and  Weapon  Delivery  Division,  Air  Force  Avionics 
Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio. 

This  report  was  submitted  by  the  author  in  June  1980. 
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SUMMARY 


The  purpose  of  this  study  was  to  generate  a  model  for  use  in  determining  the 
feasibility  of  detecting  radar  signals  beyond  the  normal  radar  horizon.  The 
mechanisms  to  be  considered  were  tropospheric  ducting  and  earth  diffraction.  Until 
recently,  models  of  tropospheric  ducts  assumed  that  the  ducts  were  horizontally 
homogeneous,  which  led  to  significant  errors  when  compared  with  experimental  results. 
Under  this  study,  ducts  are  treated  as  laterally  nonuniform  stratifications  in  the 
lower  atmosphere. 

A  Green's  function  approach  is  used  to  derive  an  expression  for  the  field 
inside  a  laterally  inhomogeneous  duct.  The  laterally  inhomogeneous  duct  is  assumed 
to  have  a  single  step  discontinuity.  The  formulation  for  the  field,  convergence 
criteria  for  the  step  size  and  the  number  of  modes  needed  for  a  solution  are  discussed. 
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SECTION  I 


INTRODUCTION 

A  solution  to  a  very  "sticky"  problem  is  given,  using  the  Green's  function 
method.  The  extension  of  the  solution  for  propagation  in  a  uniform  medium  to  a 
medium  composed  of  steps,  to  represent  the  slow  variation  in  refractive  index  with 
distance  along  the  direction  of  propagation,  allows  a  solution  for  propagation  in  a 
laterally  inhomogeneous  medium.  The  coupling  between  normal  modes  in  each  region  is 
easily  separated  out  of  the  solution  in  the  present  formulation.  The  theory  could  be 
used  to  study  the  problem  of  propagation  of  underwater  sound  waves  in  shallow  water 

with  slowly  varying  depth. ^  The  problem  of  propagation  in  a  laterally  inhomo- 

(21 

geneous  duct  was  investigated  by  Bahar'  using  an  iterative  solution  to  Maxwell's 
equations  directly. 
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SECTION  II 


ANALYSIS 

The  geometry  of  the  propagation  problem  is  shown  in  Figure  1.  The  term  duct 
refers  to  the  concept  of  the  trapping  of  modes  and  the  resulting  propagation  over 
long  distances. 


1 1 

Figure  1.  Geometry  for  single  step  discontinuity  in  a  tropospheric  duct. 
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Our  approach  will  be  to  find  the  Green's  function  for  the  bounded  waveguide 
cross  section  in  Figure  1.  In  particular,  we  will  analyze  the  effect  of  the  step 
size,  (a,,-^),  and  the  refractive  index  contrast,  n2-n-|. 

The  electric  and  magnetic  fields  will  satisfy  the  two-dimensional,  time- 
harmonic,  Helmholtz  equation 


1  iL  (r  T_  3^G  ,2g 

r  3r  3rJ  r2  2  K  b 


6(r-rs)6(<fr-<|>s) 


(1) 


0  <  <j>  <2tt 

where  (r$,4>s)  and  ( r )  are  the  source  and  observation  coordinates  respectively, 
k  is  the  wave  number,  and  the  time  dependence  is  e‘1wt.  Considering  regions  (I) 
and  (II)  separately,  G  must  meet  the  periodicity  requirement,  and  we  have 


1 1_ 

r  3r 


(r  — ) 
'  3r; 


32G 

3<})2 


k2G  = 


6(r-r) 


-4>s-2mr) 


n=-co 


(2) 


We  look  for  solutions  to  (2)  with  a  singularity  at  rs  ,<J>S  on  each  "Riemann  sheet", 
n  as 


G(r,  r')  =  y^Gjr,  r^),  r^ 


(r‘  ,  <}>s  +  2nir) 


where  G  satisfies 

(X) 
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r 


a0±r 


<  00 
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The  completeness  relation  is 


-1 

2tt  i 


gr(r’rs;X)g4 


(<t>,  4>s ; x )dA 


(3) 


(4) 


where  the  contour  (counterclockwise)  in  (4)  is  selected  to  enclose  all  the  singu¬ 
larities  in  the  complex  A-plane. 
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If  we  define 

v  =  /A"  ,  lm(v)>0 

dv  =  —  dA 

V/\ 

then  the  Green's  function  g^  on  an  "infinite"  angular  transmission  line  is 


g.U,  4>s;v2)  =  — 


i  v  ( 4>-d>  -2n7T ) 


with  | g^ | -*•  o  as  |v|-*°.  Substituting  (6)  in  (4)  gives 

i  °5.  '  o  i v(4>-<*>  -2mr) 

J^)  =  2?  f  9r(r,rs;v2)e  dv 

%oo 


with  Fourier  inversion 


go 

gr(r,  rs;vZ)  =  J  G^tr.r'  *,4>,«^s+2mT) 


-iv(<t>-<(>  -2mr) 


From  (3),  gr  satisfies 


I  d  /  Ti  ,  .  c  v 

r  >  ^  ^  +  ^  ‘  ?  9 r  ” 


6(r-rs) 


With  a  perfect  absorber  at  | <p-d>s |  =tt,  Figure  1  corresponds  to  the  case  n  =  o  for 
(7)  and  (8).  Much  of  the  discussion  up  to  this  point  can  be  found  in  the  litera¬ 
ture.^ 

The  field,  E(r),  in  region  (II)  due  to  the  aperture  field  in  the  plane  (j)=<j)o  in 
Figure  (1)  is  obtained  from  Green's  Theorem  after  integrating  over  the  cylinder  1 
at  infinity  and  over  the  aperture  plane  <p=<pQ  yielding 


E(r)  =  J 


Gtt)  !!_  i 


where  by  analogy  with  the  problem  to  the  left  of  the  aperture 


6i2^  (r,r'  ;<p,<p' )  =  ~  T  (r,r '  ;v2)eiv^_4;  ^  dv, 

•^-00 


from  which  it  follows  that 


1  °r  (2),  ,  2S  i  v  (d)-<t)  ‘ ) 

=  2^r  /  9r  (r,r  ;v  )e 

•/-oo 


We  now  make  the  parabolic  wave  equation  assumption 


-  ika,  E(r') 


Subsituting  (12)  and  (13)  into  (10)  gives 


E  [r)  =  J  ^rrE{r')  ^  gf  ;vZ)eu'^'i  dv, 


-00 


and  because  we  are  interested  in  solutions  when  v  =  ka-j ,  (14)  becomes 

oo  oo 

r /  ,  kal  f  dr'  c/  i \  f  (2) /  ,  2,  iv(0-0'). 

E(r)  =  —  /  -p-  E(r')  I  g^.  '(r,r';v  )e  v  dv  . 


Now,  E(r’)  is  the  field  in  the  plane  $=<f>0  due  to  the  point  source  at  r$,  <fs  neglecting 
reflected  fields  and  is  given  by 


E(r') 


/m  ?  iv(<i> -<U 

gj,  (r1  ,rs;v  )e  dv  . 
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Substituting  (16)  into  (15)  gives  the  solution 


oo  co  oo  •  /  \ 

E(r)  =  5/  ^  /  dV2  9‘1)(I-.rs;.f)g<2>(r.r--.vl)e1Vl<4>0'*s)- 


-00  -CO 


iv2U-<i>0) 


which  is  easily  generalized  to  2-step  discontinuities  as 


^7  ^  7 fci  / dU2  /■ 


•gi1  7r' ,rs;v2)  gj.2^(r'  ,r",v2)  gj3^ (r" ,r;v2) 


2;  yr 


i  V1  ( 4>1 )  i  ^2  ( Ci,2-<^1  ^  ^v3^2"^o^ 


The  Green's  function  for  the  source  rg  in  the  duct  is  given  by  the  "broken" 
function 

T(v)[H^}(krs)  +  R-]  (v)H^(krs)J  H<1}(kr)  ,  ai  <  r  <  - 
a(v)[H^(kr)  +  R2(v)H^2^  (kr)]* 

g(r,rs;v2)  =  -[H^1  )(krs)  +  R]  (v)  H^2)(k>"s)]  ,  rg  <  r  <  ^  (19) 


a(v)[H(])(krs)  +  Rs(v)  H^2)(krs)]- 
*[H^(kr)  +  R-|  (v)  H^(kr)] 


•  ao  1  r  1  rs 


In  (19)  and  are  Hankel  functions  of  the  first  and  second  kind. 

v  v  v 
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where  a-|  is  the  height  where  the  refractive  index  changes  from  k  to  k  where 

k  =  k(l  -  An)  | 

and  aQ  is  the  radius  of  the  earth  and  An  is  the  refractive  index  contrast. 

The  Green's  function  if  the  source  is  outside  the  duct  is 

Ia(v)[H^(krs)  +  R2(v)H^2^  (kr$)]H^ (kr) ,  r$  <  r  <  « 

a(v)[H^(kr)  +  R2(v)  H^2^{kr)]H^(krs),  a2  £  r  <_  rg  ( 

T(v)[H^(kr)  +  R-|  (v)  H^2y  (kr)]H^(krs) ,  aQ  <  r  <  a2 

where  the  height  of  the  duct  is  now  labeled  a2,  because,  as  we  see  in  Figure  2, 
the  integration  point  for  Region  II  eventually  lies  outside  the  duct. 

From  the  jump  condition 


l£L  ia  -1 

3r  3r  2tt  r 

r=r  +c  r=r  -e 


and  the  use  of  the  following  Wronskian 


H^'(kr  )H^(kr  )  -  H(1^(kr  )H^'(kr  )  =  -[p-  , 

v  v  s  v  s'  v  v  s'  v  v  s  nkr  ’ 

s 


we  obtain 


a(v)  = 


8[R2(v)-R1(v)] 


and  the  "resonance"  condition;  i.e.,  R2(v)  =  R^ (  v).  From  the  boundary  conditions 
at  r=a^  ,  i.e.. 


9(r,rs;v  ) 


=  g(r,rs;v  ) 


r=a-|-£ 


1  r=a-j+c 


!r=a^-e  '^a^e 
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we  find 
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R2(v)  =  -{k  H^),(ka2)[H^1)(ka2)  +  R]  (v)H^2)  (ka2)] 

-  k  hJ1}  (ka2)[H{1)'(ka?)  +  R. (v)H(2) ' (ka2)]}  / 

v  2  v  2  1  v  2  (34) 

{k  H<2)l(ka2)[H^1)(ka2)  +  R^v)  H<2)  (ka2)] 

-  k  H<2)  (ka2)[H(1)l(ka2)  +  R1(v)H^2)l(ka2}]} 

and 

T(v)  =  a(v)[H(1)(ka2)  +  R2(v)H{2)(ka?)]  / 

v  2  2  v  2  (35) 

[H^1  )(ka2)  +  R1(v)H|)2)(ka2)]. 

We  now  turn  our  concentration  to  the  aperture  integration  aQ  _<  r'  <“  in  (17).  The 
geometry  of  the  problem  is  given  in  Figure  2,  where  we  will  confine  our  analysis  to 
the  case  where  aQ  <_  rg  £  a-|  and  aQ  <  r  <  a?. 


r 


Figure  2.  Geometry  for  aperture  integration.  Modes  in  regions  (I)  and  (II) 
are  and  v2,  respectively. 


> 
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Using  the  "broken"  functions  in  (19)  and  (20)  with  appropriate  interpretations  for 
the  source  and  observations  points  (i.e.,  in  the  aperture  plane  the  integration 
point,  r‘ ,  becomes  the  observation  point,  r,  in  Region  (I)  while  the  integration 
point,  r',  becomes  the  source  point,  rs,  for  Region  II),  we  have 


oo  r 

/^9^(r,,rs;v?)g^(r,r,^)  =  a(v1  )ot(v2)<j>  (krs)^(kr)  f  (kr' )6v^(kr‘ 
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where 


$v(kr)  =  } (kr)  +  R] (v)H<2) (kr) 

ipv  ( k  r )  =  H^1  ^  (kr)  +  R2(v)H<2)  (kr) . 


(37) 


From  the  differential  equation  (9)  for  the  radial  functions  and  linear  combinations 
as  given  in  (37)  we  have 


f  ^  Vkr,)Vkr,)  =  7TL?7cS(kr) 

“a  2  'W  £ 


'l 


3r’ 


B<p, 


V  r 


r  =r 


‘\(kr)  IF 


a* 


T~T7  [Kjka o)  ~aF 


(v{-v|)  2 


a<{> 


] 


(38) 


-  K  (kaft) 


r'sa_ 


^'^o'  aF 


r’*a. 


a* 


36 


FT\  ~Brr 


(vf-v‘)  v2 


r’*r 


(kr) 


r  *r 


10 


because  *  (kr)  satisfies  (28).  Also 
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2  27  [  V(krs)  IF 


(vj-\£)  2 


9^, 


-<K,  (krc 


Vo 


v, v  s'  9r 

r’  =  rs 


]. 


The  following  integral  involves  both  the  wave  numbers  k  and  k;  i.e., 


iO) 


a  8h(1v  ^kr,) 

'a,  •  HV(kr,)%{kr,)  =  7T-T:  t%(ka2}  — §F 

1  2  (vrv2)  2 


fd-f 


dip 


"Hv,  ka2  ~Jrr 
r’=a2  1 


’1 


,  2  2i  fV(kal)  9r 1 

( v-|  ”V2 )  2 


9H(1^  (kr*) 
V1 


r  -a. 


9tjj 


]  + 


2Ank 


,  2  2% 
r'-a,  (vl'v2} 


] 


2 

(41) 


/ 


d'dr*  H^(kr'  H  (kr1) 
V1  v2 
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and,  finally. 


r  a  3H\^(kr') 

-  -  -rV  [  — L, -  -H<’>(ia2) 

*  lW  1  r'=a9  1 

aH^^kr') 

-  ]  • 


Substituting  (38),  (39),  (40),  (41)  and  (42)  into  (36)  gives 


oo 


gj1)(r,,rs;v12)gj2)(r,r';v|)  =  (  p-  ) 


[R2(v2)-R1  (v2)][R2(v1  )-R1  (v^Kvj-Vj) 


3<t> 

'>v1(krs)‘tv](k'")r  tV(kr)  ~~j£  -%<kr)  -IF-  ] 

r'=r  r ' =r 


30  3ty 

+  %(krK2(krs)rs  [  '"v1(krs)  IF  -%(krs}  IF  3 

r’=rs  r1 =r$ 


♦  t*Vi(krs)  -jp-  -  *v,<kal> 


v ,(krs> 


3H(1 ^(kr‘ ) 
V1 


r,=a1  V1  '  ^kal> 


V(kal) 


3H(1)(kr’), 


*  *v2(kr>  *v ,(krs>  H-(TT,T— '  >2  K  <k»2>  — 5F -  -Hi!,<ka2>  1 

v-|  '  r  r'=a2  r'=a2 
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(5) 


K  (kai^  2 

+0  (kr)<f>  (kr,)  — ttt-  — -  2Ank2  r'dr'  H^(kr'H  (kr1) 

V  V  s  H^tkaJ  a  V1  v2 

V1  1  al 


^  ( ka-j )  [^(vgJ-R, (vgJlCH^1  ^(kagjH^1  ^  ‘  (kagJ-H^1  } ' (ka2)3 

+4  (kr)*  (kr  )  ,]  > . a .  (4i/ir)  - - - - - - - - - - - j 

2  1  %  (kal}  [k  H^)‘(ka2)<i)V2(ka2)-hH^)(ka2)0>vMka2)]  I 


We  will  show  that  terms  (1)  and  (2)  in  (43)  represent  the  "uncoupled"  or  zeroth  order 
solution  and  (6)  represents  the  first  order  coupling,  while  the  remaining  terms; 
i.e.,  (3),  (4),  and  (6)  are  all  of  second  order.  Consider  the  following 
combination  of  terms  from  (3)  and  (5);  i.e., 


(ka-j)  3H^(kr') 

^(kr)%(krs)  )  [a2%(ka2}  “IF1 - 

V1  1 


3Hj1}(kr') 

‘al'kv2(kal)  ~SF -  _  ] 


V(kai} 

=  k<*>  (kr)*  (kr  )  n  . - 

v2  V1  s  H'  (ka, ) 
V1 


[Vv2(ka2)H!j)  (ka2)  ■  al^v2(kal)H!11)  (kal)]- 


We  will  make  use  of  the  following  asymptotics 

m  2  1/3 

(vz)  ~  (  ^  )  Wi  (t  +  x) 


(?)  ?  1/3  (2> 
H(2,(vz)  ~  (  )  Wi  (t  +  x) 


where  Wi^  and  Wi^  are  Airy  functions  with  ^ 

t  -  (^)2/3[l-  (^M  or  (v  =  ka  +  (  IjS.  )V3  t) 


and  the  dimensionless  "radial  wavenumber" 


x  =  ka/ (  *=■  ) 
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Using  (45),  (46)  and  (47)  we  find  (44)  becomes 

,  o  V3  0)  (2) 


1/3  (1) 


k  [Wi  (t2+Xr)+R1(t2)Wi  (t2+Xr)](^-)  [Wi  (t1+Xs)+R2(t1)Wi  (t1+X$)]. 


1/3  (1) 


?  1/3  (1) 

^  \  Hi 


a2  (^~)  [Wi  (t2)  +  R2(t2)Wi  (t2)]  Wi  (t]  +  xa  +  xQ) 

2  ka0 


P  1/3  (1)  (2)  ?  1/3  (1 ) 1 

'al  [Wi  (t2+Xa)+R2(t2)Wi  (t2+xa)]  (4~)  Wi  (t 

1  ka  -j 


1  +  xD)j 


ka  1/3 
I  \ 


xr  =  k  ( a i -r ) / (— ^ — )  > 


ka,  1/3 

xfl  =  k  (a2-a i )/ (  jjr  )  , 

ka,  1/3 

xs  =  k ( a -j  -  r$)/(-y-) 


ka,  2/3 

Xp  =  2An(-^— ) 


Now,  we  see  if  the  radial  wavenumber,  xa,  is  small,  (48)  becomes 

0 

„  , p  1/3  p  1/3  .  1/3  (1)  (2) 

k(aJ/3-a]/3)(y  (^f-)  (4)  [Wi  (t2+xr)+R1(t2)Wi  (t2+xr)]. 


•[Wi1)(t1+xs)+R2(t1)Wi2){tl+xs)]  CWl  (t2)+R2(t2)Wl  (t2)]  * 

(1)’ 

•WI  (t]  +  x0) 


which  Is  negligible  compared  with  (1),  (2)  and  (6). 
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Also  from  (3)  and  (5),  we  find 


%(kr)V,(krs)al  [%(kal}  TF 


-%(ka1}  ~W 


o  2/3  (1)  (2) 

-  2ikal(|)v^(kr)0v^  (krs)R2(t2)  (^-)  [Wi  (t]  )Wi  (t2  +  xg)  -  (50) 

(2)  (1) 

-  Wi  (t-j  )Wi  (t2  +  xg)] 

which  is  negligible,  again  compared  with  (1),  (2)  and  (6).  Consider  (6),  and 
the  rational  function 

- } - - - £ -  .  (51) 

[k  (ka2)^(ka2)-  k  (ka2)<0^  (ka2)] 

(1)  (1)'  (1)  (1)’ 

[Wi  (t2+xD)Wi  (t]+xD+xa)-Wi  (t1+xD+xg)Wi  (t2+xD)] 

:  UF  ro  m  rn  nr^  jtt 

k  Wi  (t2+xD)[Wi  (t2)+R1  (t2)Wi  (t2)]-k  Wi  (t2+xD)[Wi  ( t2 )+R-j  ( t2 )Wi  (t2)] 

which  is  negligible  for  small  radial  wavenumber,  xg.  Note,  even  if  x^  =  0,  the 

(1)'  (2) 

denominator  in  (52)  is  non-zero  because  it  involves  the  Wronskain  Wi  (t9)Wi  (t9)- 
(1)  (2) 1  L  L 
Wi  (t2)Wi  (t2)=-2i/iT.  Of  the  terms  in  (43)  remaining,  (1)  and  (2)  represent  the 

uncoupled  normal  modes  and  are  given  by 

9  4>v 

^v1(krs)V(kr)r  [0v9(kr)  IF  '  \(kr)  “§F  ] 

i  ’  L  r'=r  L  r'=r 


(4i/7i)[R1  (v2)  -  R2(v2)]  ^(krs)  0y^(kr) 
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d<t> 

K2(kr,K2<krs,rs  [*v,(krs>  TF 


dip, 


-  <*>,  (krj 


1 


v,v  s'  9r 

r'=r 

s 


r  =r_ 


(53) 


=  ( 4 i / tt )  [R2(v])  -  R-j  (v-j )]  <t>^(kr)^(krs) 


where  use  of  (23)  is  made.  The  term  involving  the  coupling  between  regions  (I)  and 
(II)  is  (5)  is  given  by 

a. 


(ka-j )  °2 

( kr ) d>  (krc)  — - -  2Ank^  f  r'dr'  H^^(kr')^  ( kr ' ) 

/  '1 
dl 


V-|  1 


(54) 


Returning  to  (17),  we  have  the  v-|-  and  v2~integrations  remaining.  These  are  easily 
performed  using  Cauchy's  Theorem;  the  integrands  are  analytic  except  at  the  simple 
poles  which  are  solutions  of 


v,  v, 

l'v2  2'v2 

The  residues  at  the  simple  poles  are 

(1) 


xD[Wi  (t  +  xn)]‘ 


3](t)  =  (tt/2i  )  |  ("2")  rrr  rzr  rn 


[Wi  (t)Wi  (t+xn)-Wi  (t)Wi  (t+xn)]  [Wi 


!2Wn2 


-1 


(55) 


(56) 


with 


xQ  =  k ( a-j  -  afl)/(ka1/2) 
2  2 


1/3 


(57) 
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The  v-|-  and  v2-integrations  yield  the  desired  result 


E(r)  =  (-1/32)  - Wt,-  Va^a^kr^fkr) 

(ka/2)l/J  ^1  1 

•exp[i ) ( ka-j  +  i  (ka^'2)1/3  t™  ] 

+ - y~'ai(t2)<}>  m(kr)\fj  m(kr  ) 

(ka/2)l/J  ^  1  1  h  *2  s 


with 

2  2  (  ka  a  ) 

[(v*)  (vj)  ]  »  2k3]  k(aia2)  +  (-/)  [tf(^)  t$]j  (59) 

In  (58),  E(r)  is  the  field  normalized  by  the  source  intensity,  this  result  is  used 
for  numerical  experimentation  and  has  been  further  normalized  as  W  -  |E(r)|  (-1/32). 
Therefore,  W  is  dimensionless. 
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SECTION  III 


REMARKS 

1)  The  single  sums  over  the  modes  in  (58)  represent  the  field,  assuming 
uncoupled  normal  modes^’^.  The  assumption  proceeds  by  taking  the  boundary 
conditions  to  be  independent  of  the  coordinate  in  the  direction  of  propagation,  but 
the  boundary  conditions  in  the  normal  direction  are  the  same  that  would  be  applied 
for  perfectly  stratified  media.  It  is  also  assumed  that  eigenfunctions  correspond¬ 
ing  to  a  particular  normal  mode  are  orthornormal ;  i.e., 

oo 

/  ^-t  (kr')d>  (kr‘)  =  6 

J  r  nv  /rm  '  mn 

ao 

Unfortunately,  in  many  theories,  it  is  often  very  difficult  to  justify  when  to 
neglect  the  coupling.  The  solution  given  in  (58)  allows  a  direct  determination  of 
the  effect  of  coupling.  In  (58)  if  the  coupling  is  negligible,  the  solution  suggests 
the  angular  position  of  the  step  is  unimportant,  and  one  must  only  remember  which 
medium  he  is  in;  e.g.,  source  in  medium  (I),  observer  in  medium  (II). 

2)  Cho  and  Wait^  gave  a  derivation  for  the  fields  in  a  stepped  model  for 

a  non-uniform  duct  which  employed  the  use  of  a  non-Hilbert  space  inner  product; 
i.e.,  <  <J>n ,  <|>  >  ,  instead  of  the  usual  definition  in  terms  of  a  complex-valued 

function  or  ordered  pairs  with  inner  product  <<)>,<(>*>.  The  natural  metric 

r  Tn  m 

1/2 

{x-y,  x-y} 

is  a  real  nonnegative  quantity  and  represents  the  physical  quantity  power.  Recall¬ 
ing  a  metric  space  is  complete  if  every  Cauchy  sequence  is  a  convergent  sequence, 
the  usual  definition  of  a  Hilbert  space  is  an  inner  product  space  which  is  complete 
with  respect  to  its  natural  metric.  The  Cho  and  Wait  result  can  be  explained  by 
the  use  of  "biorthogonal "  coordinates^.  Let  (vn)  be  the  set  of  nonzero  eigen¬ 
values  of  the  differential  operator 

£  _i  (  r  —  )  +  rk2 

and  let  {<t>n >  be  the  corresponding  eigenfunctions.  The  nonzero  eigenvalues  of  the 
adjoint  (formal)  operation  £  are  given  by  {v*}  and  the  corresponding  eigenfunctions 
will  be  denoted  by  Now,  take  for  the  set 

{>km>  =  {<(>!  <P2  »•  *  •) 
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Then,  indeed,  the  innner  product  will  satisfy 


In  fact,  Cho  and  Wait's  result  for  <  >  equals  our  result  in  (56);  i.e., 

<  ‘‘V^m  >  =  ai^^-  This  yields  the  interesting  conclusion  that  the  Cho- Wait  inner 
product  will  equal  zero  if  and  only  if  a-|(t)  equals  zero  which  requires  the 
existence  of  a  double  root!  This  can  occur  even  in  a  single  section  duct.  Proof: 
Since  the  denominator  of  a^(t)  is  a  rational  function,  the  only  singularities  it 
can  have  in  the  entire  complex  tplane  are  poles.  A  double  root  suggests  degenerate 
modes  in  the  two  regions;  i.e.,  the  spatial  distribution  of  sources  across  the 
aperture  plane  <J>=4>0  has  the  same  wavelength  as  the  normal  modes  being  driven  and 
a  resonance  occurs.  The  integral  formulation  used  here  in  terms  of  the  Green's 
function  approach  completely  sidesteps  the  issue  of  what  the  normalized  "eigen¬ 
functions"  should  be  in  the  biorthonormal  coordinate  approach.  It  may  be  that  since 
the  residue  in  the  Green's  function  approach  equals  the  inner  product  in  the  bio- 
rthornormal  case  and  because  the  H^(kr)  are  dense  in  our  Hilbert  space,  any  function 

Vm  y*  (i \ 

can  be  approximated  to  within  e>o  by  ^-<am  '(kr).  The  problem  comes  in  finding  out 

m 

how  to  express  the  am's.  This  is  probably  an  example  of  a  problem  where  the  solution 
can  be  found  by  a  Green's  function  method  but  only  a  generalization  of  the  notion  of 
eigenfunctions  permit  a  solution  in  terms  of  the  latter.  The  other  point  is  that  the 
residues  come  out  naturally  in  the  Green's  function  method. 

3)  The  double  sum  in  (58)  depends  upon  the  location  of  the  vertical  step, 

<t>  ,  and  represents  the  coupling  from  mode  m  to  mode  r..  The  magnitude  of  this  term 
depends,  on  the  electrical  step  size,  k^-a^),  as  seen  in  (58)  and  (59). 
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SECTION  IV 


EXAMPLE 

The  numerical  results  for  the  four  "height-gain"  curves  in  Figure  3  correspond 
to  the  following  choice  of  parameters: 

aQ  =  6378  km 
a-|  =  6379  km 

a2  =  6379000,  6379010,  6379020,  6379030  m 

rs  =  6379  km 

f  =  100  MHz 
An  =  25  N-units 

a0(V4>s)  =  100  km 

a0(4ls_<*)  )  =  100  km 

6  =  0.3  -  i 4x1 0- ^  (a  =  .001  Seimens/m,  er  =  10) 

From  Figure  3  it  appears  that  a  step  size  of  about  20m  causes  significant  change  in 
the  height-gain  pattern.  We  will  refer  to  this  as  the  "resonant"  step-size.  This 
would  correspond  to  a  radial  wave  number  from  (47)  of  about  0.463  radians  ( i . e . , 
about  t/8).  The  second  limiting  criteria  for  our  solution  in  (58)  is  the  number  of 
modes  required  for  convergence  of  the  series.  For  the  example  in  Figure  3,  10  modes 
gave  two  significant  figures.  The  convergence  of  the  series  is  dominated  by  the 
exponential  terms  in  the  series  for  small  m  and  by  the  asymptotic  decay  of  the 
residues  for  large  m;  i.e., 

a, (t  )  ~  exp(-  4/3  r!/2)/4  tl/Z 
I  m  r '  m  '  m 

where  for  tm  <  xQ,  the  imaginary  part  of  tm  becomes  small.  In  Figure  4,  the  effect 
of  the  number  modes  is  shown  for  a  10  km  separation  between  source  and  step  and  step 
and  observer.  At  this  distance  30  modes  are  required  for  convergence. 

In  Figures  5  and  6,  thechoice  of  the  parameters  is  the  same  as  in  Figure  3 
except  the  frequency  is  300  and  60  MHz  respectively.  Figure  7  is  the  same  as 
Figure  4  except  the  refractive  index  contrast  is  50  N-units.  The  limiting  step  size 
for  this  case  is  about  10  m.  The  resonant  step  size  for  300  MHz  is  about  5  m  for 
An  =  25.  In  Figures  5,  6  and  7,  10  modes  provided  adequate  convergence  of  the  sums. 


20 


of  Modes 


Figure  4.  Convergence  of  solution  with  increasing  number  of  modes.  Distance 

between  source  and  step  and  step  and  observer  is  10  km.  The  frequency 
is  100  MHz  and  the  refractive  index  contrast  25  N-units. 


Height-gain  run  for  300  MHz.  W  is  the  normalized  signal  strength 
The  refractive  index  contrast  is  25  N-units.  100  km  separation 
between  source  and  aperture  plane. 


(W)  IM&!9H 
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Figure  6.  Height-gain  run  for  60  MHz.  W  is  the  normalized  signal  st 
The  refractive  index  contrast  is  25  N-units.  100  m  separa 


Figure 


SECTION  V 


CONCLUSIONS 

A  Green's  function  approach  is  used  to  examine  the  effect  of  varying  the  step- 
height  in  a  tropospheric  duct  with  a  single  step  discontinuity.  If  the  electrical 
height  of  the  step  is  less  than  the  "radial"  separation, 

k(a2  -  )/(ka/2)1/3 

the  vertical  distribution  of  the  field  strength  agrees  with  the  fields  in  a  duct 
with  no  discontinuity.  This  agrees  with  a  result  obtained  by  Wait  and  Spies  for 

/  Q  \ 

an  ionospheric  ductv  .  The  location  of  the  step  in  relation  to  the  source  and 
observer  determines  the  number  of  modes  required  for  convergence. 


SECTION  VI 
REFERENCES 


(1)  A.  Nagl,  H.  Uberall,  6.  L.  Zarur,  and  A.  J.  Haug,  Adiabatic  mode  theory  of 
underwater  sound  propagation  in  a  range  dependent  environment,  J.A.S.A., 

Vol .  63,  p.  739,  (1978). 

(2)  E.  Bahar,  Propagation  of  VLF  radio  waves  in  a  model  earth-ionosphere  waveguide 
of  arbitrary  height  and  finite  surface  impedance  boundary:  Theory  and  experi¬ 
ment,  Radio  Science,  1,  (New  Series),  Vol.  8,  pp.  925-938,  (1966). 

(3)  L.  B.  Felsen  and  N.  Marcuvitz,  Radiation  and  Scattering  of  Waves,  Prentice-Hall, 
Inc.,  Englewood  Cliffs,  New  Jersey,  pp.  306-314,  (1973). 

(4)  R.  H.  Ott,  Roots  to  the  model  equation  for  EM  wave  propagation  in  a  tropospheric 
duct.  Jour.  math.  Physics,  Vol.  21,  pp.  1256-1266,  No.  5,  May,  (1980). 

(5)  J.  B.  Keller,  Progress  and  prospects  in  the  theory  of  linear  wave  propagation, 
SIAM  Review,  April,  Vol.  21,  No.  2,  pp.  234-235  (discussion  of  "Propagation 
and  normal  modees"),  (1979). 

(6)  Allan  D.  Pierce,  Extension  of  the  method  of  normal  modes  to  sound  propagation 
in  an  almost-stratified  medium.  Jour,  of  the  Acoustical  Society  of  America, 

Vol.  37,  No.  1,  pp.  19-27,  Jan.,  (1965). 

(7)  S.  H.  Cho  and  J.  R.  Wait,  EM  wave  propagation  in  a  laterally  non-uniform 
troposphere,  EM  Report  #1,  CIRES,  Boulder,  Colorado,  June,  (1977).  A  summarized 
version  is  given  in  S.H.  Cho  and  J.R.  Wait,  Analysis  of  Microwave  Ducting  in 

an  Inhomogeneous  Troposphere,  Pure  and  Applied  Geophysics,  Vol.  116,  pp.  1118- 
1142,  (1978). 

(8)  I.  Stakgold,  Boundary  Value  Problems  of  Mathematical  Physics,  Vol.  1,  The 
MacMillan  Company,  London,  p.  257,  (1964). 

(9)  J.  R.  Wait  and  K.  P.  Spies,  On  the  calculation  of  mode  conversion  at  a  graded 
height  change  in  the  earth-ionosphere  waveguide  at  VLF,  Radio  Science, 

Vol.  3  (New  Series),  No.  8,  pp.  787-791,  August,  (1968). 


27 

*U.S.Oov«rnm«nt  Printing  Offlcti  19(0  -  *87-0(4/134 


